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ABSTRACT 

We present a comprehensive study of the star-formation properties of the infrared luminous galaxy 
Arp 299 (IC 694 + NGC 3690). The observations include HST/mCMOS imaging and MMT optical 
and near-infrared spectroscopy together with HST archival data. We correct the galaxy parameters for 
extinction and use the results as boundary conditions for evolutionary starburst models. 

These models and other arguments show that Arp 299 has been going through a broad variety of 
interaction-induced star formation for the last ^ 15 Myr. In addition to nuclear starbursts in the two 
colliding galaxies, there are nearby luminous star forming regions that may be the result of molecular 
clouds breaking up and starting to form massive stars as they approach the nuclear potential. Two 
regions near the interaction region have very recently formed massive stars 4 Myr old). One of these 
regions in particular appears to have undergone a starburst of very short duration, in contrast with the 
5-10 Myr durations that are typical of nuclear events. These regions will probably form gravitationally 
bound new dwarf galaxies, although they are close enough to the more massive original galaxies that 
they will probably eventually be subsumed into them. In addition, we find 19 H Ii regions at least as 
energetic as 30 Doradus, and 21 older star clusters that likely are the products of similar Hii regions 
formed in the past. This population of supergiant Hii regions is unprecedented in normal galaxies and 
emphasizes that the effects of the interaction propagate through the entire pair of galaxies. 

Subject headings: Galaxies: individual (Arp 299, IC 694, NGC 3690) - galaxies: nuclei - galaxies: 
photometry - galaxies: stellar content - galaxies: interacting - infrared: galaxies 



1. INTRODUCTION 

Ultraluminous infrared galaxies that rival quasars in en- 
ergy output were first identified nearly thirty years ago 
(Rieke & Low 1972), and were shown by ZR^5 fifteen years 
ago to exist in numbers comparable to those of quasars in 
the nearby Universe (e.g., Sanders et al. 1988). The exact 
nature of these extreme objects has been of great inter- 
est and debate since. Among the ultraluminous infrared 
class (LiR > 10^^ L0) it is found that a large percentage 
are interacting/merging systems containing active galactic 
nuclei. The process of merging with accompanying super- 
starbursts appears to be an important stage in galaxy evo- 
lution, possibly even converting spiral galaxies into ellip- 
ticals. Sanders et al. (1988) suggested that the infrared 
luminous phase is the initial stage for the appearance of a 
quasar (see Sanders & Mirabel 1996 for a recent review). 

To understand these objects requires observations at 
high physical resolution. However, because their space 



density is low at the current epoch, there are few nearby 
examples where such resolution is possible. Arp 299 
(NGC 3690/IC 694 or Mrk 171) is of high infrared lumi- 
nosity (LiR = 5 X 10^^ Lq), placing it among the luminous 
infrared galaxies, but also is one of the nearest examples 
of interacting starburst galaxies (distance D = A2 Mpc for 
Ho = 75km s'^ Mpc"!). 

Within Arp 299, a number of bright infrared and radio 
sources appear to be the regions of star formation and/or 
nuclear activity responsible for the large overall luminosity 
(see Gehrz, Sramek, & Weedman 1983; Wynn-WiUiams 
et al. 1991 and references therein)^. As in other lumi- 
nous and ultraluminous galaxies, there is a high concentra- 
tion of molecular hydrogen within relatively small regions 
around these sources. From CO maps, Sargent & Scoville 
(1991) estimated densities of molecular gas ~ 8 x 10^ Mq 
pc^'^ in IC 694, ~ 3 x 10'* Mq pc~^ in components Bl and 
B2 of NGC 3690 and ~ 2 x lO'' Mq pc'^ at the interface 
of both galaxies (region C-l-C'). These general characteris- 



^We use the notation introduced by Gehrz et al. (1983) for the interacting pair of galaxies Arp 299: the nucleus of IC 694 (eastern 
component) is source A, and the sources in NGC 3690 (western component) are Bl, B2, C and C' (see also Wynn-WiUiams et al. 1991 and 
Alonso-Herrero, Rieke, & Rieke 1998). 
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2 The interacting galaxy Arp 299 

Table 1 

Log of the NICMOS observations 



Camera Filter tgxp Field of view 



NIC3 


F164N 


1024 




51.2" 


X 51.2" (IC 694- 


+NGC 3690) 


NIC3 


F166N 


1024 




51.2" 


X 51.2" (IC 694- 


+NGC 3690) 


NIC2 


F160W 


208 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F222M 


528 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F237M 


864 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F187N 


320 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F190N 


320 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F212N 


2560 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NIC2 


F215N 


2560 


19.5" 


X 19.5" 


(IC 694), 19.5" 


X 19.5" (NGC 3690) 


NICl 


FllOM 


256 


42.6" X 21.0" 1 


(IC 694 + NGC 3690) 8 image mosaic 



tics are consistent with numerical simulations of collisions 
between gas-rich galaxies (see Barnes & Hernquist 1996 
and references therein) that show how interactions lead to 
transportation of large quantities of molecular gas into the 
centers of galaxies. This concentration of gas leads to a 
strong burst of star-formation and may activate an AGN. 

However, these properties make the centers of activity 
in luminous and ultraluminous infrared galaxies virtually 
impenetrable in the visible and ultraviolet. Wc report near 
infrared HST/NICMOS images that penetrate much of the 
interstellar extinction and arc at an unprecedented level 
of resolution for this spectral region, up to 0.1" corre- 
sponding to 20 pc. We also use WFPC2 images of the 
system together with ground-based optical and infrared 
spectroscopy to derive the star formation history. These 
new data are placed in the context of the characteristics 
of the galaxy measured in other spectral regions, as sum- 
marized briefly in the preceding paragraph. 

2. OBSERVATIONS 
2.1. NICMOS observations 

HST/mCMOS observations of Arp 299 were obtained 
on November 4 1997 using all three cameras. In Table 1 
we give the log of observations, where column (2) is the 
filter, column (3) is the integration time and column (4) 
is the field of view of the images. The observations were 
obtained in a spiral dither with a 5.5 pixel spacing and 
with two, three or four positions. The plate scales for 
NICl, NIC2 and NIC3 are 0.045" pixel"!, 0.076" pixel"! 
and 0.20" pixel"! respectively. 

The images were reduced with routines from the pack- 
age NicRed (McLeod 1997). The main steps in the data 
reduction involve subtraction of the first readout, dark cur- 
rent subtraction on a readout-by-readout basis, correction 
for linearity and cosmic ray rejection (using fullfit), and 
flatfielding. Darks with sample sequences and exposure 
times corresponding to those of our observations were ob- 
tained from other programs close in time to ours. Usually 
between 10 and 20 darks were averaged together (after 
the subtraction of the first readout) for a given sample 
sequence. However, in the case of FllOM we got better 
removal of the artifacts using a dark-|-sky image generated 
by median combining all the dithered galaxy images for 
this filter. Flatfield images were constructed from on-orbit 
data. Since our images were obtained after August 1997, 



no correction for the pedestal effect was necessary. For 
more specific details on the reduction of the NIC3 images 
see Alonso-Herrero et al. (1998). 

After removal of artifacts, the dithered galaxy images 
were registered to a common position using fractional pixel 
offsets and cubic spline interpolation, and combined to 
produce the final images. The NICl FllOM images were 
mosaiced. Since there are no bright common sources in the 
overlapping regions of the images, we used the world co- 
ordinate information contained in the header to compute 
the offsets. All the NICMOS images were rotated to the 
usual orientation north up, east to the left. 

Prior to flux-calibrating the images, we measured the 
background on blank corners of the images. Although 
this interacting pair of galaxies is known to be extended 
beyond the field of the NIC2 images, the background 
values mcasrued in the NIC2 F222M fAtci image are in 
excellent agreement with the background measurements 
taken during the Servicing Mission Observatory Verifi- 
cation (SMOV) program. The fiux calibration was per- 
formed using the conversion factors based on measure- 
ments of the standard star P330-E (Rieke 1999). The 
fully-reduced narrow-band images (NIC2 F187N, F190N, 
F212N and F215N, and NIC3 F164N and F166N) were 
first background subtracted, then fiux calibrated, and fi- 
nally shifted to a common position. The adjacent contin- 
uum image was subtracted from the corresponding line + 
continuum image to produce the final continuum-free line 
images. 

Figure 1 shows the NIC3 F166N image of the pair of in- 
teracting galaxies. The NIC3 F166N filter at the redshift 
of Arp 299 contains the emission line [Fell] 1.644 /im. The 
total field of view is 51.2" x 51.2". The corresponding im- 
ages of the continuum subtracted [Fe ll] 1.644 /xm line emis- 
sion of the system can be found in Alonso-Herrero et al. 
(1998). Figure 2a and Figure 2b show both components 
through different filters as indicated in the figure caption. 
The field of view of all these images (except for the NICl 
FllOM inset in Figure 2a) is 19.5" x 19.5". The NIC2 
F222M and F237M images are similar in morphology and 
are not shown here. 

2.2. Optical and ultraviolet images from the HST archive 

A Wide Field Planetary Camera 2 (WFPC2) image of 
Arp 299 was retrieved from the HST data archive. This 
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Fig. 1. — NIC3 F166N image ([Fc ii] 1.644 /im + continuum) of the interacting system Arp 299 on a logarithmic scale. The total field of 
view is 51.2" x 51.2". Orientation is north up, east to the left. 



Table 2 

Low-resolution /C-band infrared spectroscopy for a 1.2" x 1.2' 

EXTRACTION APERTURE. 



Source 


/(H2) 


EW(H2) 


/(Br7) 


EW(Br7) 


(CO)sp 


(CO)ph 




(erg cm ^ s ^) 


(A) 


(erg cm ^ s ^) 


(A) 






Bl 


9.52 X 10*1'^ 


10.2 


8.56 ±0.5 X 10- 


9.0 ±0.5 


0.14 


0.07 


B2 


< 2.4 X 10-15 


< 4 


< 3.1 X 10-15 


< 5 


0.31 


0.17 


C 


2.95 X 10-15 


5.1 


2.43 ±0.1 X 10-" 


41 ±2 
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Fig. 2a.— IC 694. Top panels, left: NIC2 F160W (//-band continuum), and right: NICl FllOM (continuum at 1.1 /im) right. Middle pan- 
els, left: optical WFPC2 F606W, and right: continuum-subtracted NIC3 F166N ([Fe n]1.644/im). Bottom panels, left: continuum subtracted 
NIC2 F215N (H2 at 2.12 ^m), and right: continuum-subtracted NIC2 F190N (Paa). Orientation is north up, east to the left. The field of 
view of all the images except the inset is 19.5" x 19.5". 




Fig. 2b.— Same as Figure 2a but for NGC 3690. 
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image was taken through the F606W filter on September 
17 1994 with a total integration time of 500 s and spatial 
sampling of 0.046" pixel" ^. Standard pipeline reduction 
procedures were applied. The image was rotated to the 
usual orientation north up, east to the left. Since a sin- 
gle integration was taken for this galaxy, the cosmic ray 
removal is problematic. The cosmicray task in iraf suc- 
cessfully eliminated the most point-like hits but artifacts of 
extended hits remain in the data. The image is presented 
in Figure 2a and 2b, with the field of view matching that 
of NIC2. 

We also retrieved UV FOC images from the if5T archive 
for each of the components in the system, taken through 
the UV F220W filter with integration time of 897 s for each 
component of the pair. Standard pipeline reduction proce- 
dures were used. Again the images were rotated to match 
the orientation of the NICMOS images. The plate scale 
of these images is 0.022" pixer^ The image of NGC 3690 
was part of the study by Meurer et al. (1995) of the UV 
properties of starburst galaxies. 

2.3. Ground-based Ha imaging 

CCD Ha narrow-band images (on-line and off-line) of 
Arp 299 were obtained from the La Palma Observatory 
archive. The images were taken with the 1.5 m JKT tele- 
scope on January 21 1990 with a CCD-RCA2 detector. 
The integration time for the continuum-|-line and contin- 
uum images was 1,800 s each. The reduction steps involved 
bias removal, dark current subtraction, and flat-fielding 
with both dome and sky flats. The plate scale of the im- 
ages was 0.4"pixel~^, providing a total field of view of 
130" X 210". 

Since no images of suitable standard stars are provided 
by the La Palma archive, the final continuum-subtracted 
Ha image was calibrated with the spectroscopic data from 
Armus et al. (1989) and Gehrz et al. (1983). The latter 
reference provides fluxes through an effective aperture di- 
ameter of 5 ". The former authors provide the Ha fluxes 
for components A and C of the interacting pair through 
a 2" X 4.5" aperture. Because of difficulties in setting the 
slit accurately on the heavily obscured sources, there are 
significant uncertainties in the line measurements. For ex- 
ample, Armus et al. (1989) report source C to have twice 
as much flux as source A, whereas our Ha continuum- 
subtracted image indicates that source C is four times as 
bright as source A through the same aperture. Our flux 
calibration was based on an average of the values from 
the spectra; the effects of the uncertainties will be dis- 
cussed when they are relevant to the analysis. In Figure 3 
we present images of both the adjacent continuum and 
contimmm-free Ha line emission, showing the positions of 
the slits used for both the optical and near-infrared spec- 
troscopy (Sections 2.4 and 2.5). 

2.4. K-band spectroscopy 

High (A/AA ~ 3000) and intermediate (A/AA ~ 1000) 
resolution spectra centered at wavelengths 2.32 /xm and 
2.20 yum respectively were obtained at the Multiple Mir- 
ror Telescope with the FSPEC infrared spectrometer 
(Wihiams et al. 1991) on February 11 and 12 1998. The 
slit size was 1.2" x 30" with plate scale 0.4"pixer^ The 
slit was centered on Bl at position angle PA = 145° to 
cover both sources Bl and B2, and at PA= 80° centered 



on source C. The integration times were 96 minutes for 
soiirccs Bl and B2, and 28 minutes for source C for the 
high-resolution data, and 14 minutes for each of the three 
sources for the low-resolution data. 

Observations were obtained for each galaxy at three or 
four positions along the slit, integrating for 2 or 4 min- 
utes at each position. Standard stars were measured in 
a similar fashion (using shorter integrations), interspersed 
with the galaxy observations. The standards were selected 
to be at similar air masses as the galaxies so they could 
be used to correct for atmospheric absorptions. The data 
reduction process involves dark current subtraction, flat- 
fielding, and sky subtraction. The correction for atmo- 
spheric absorption is performed by dividing the galaxy 
spectrum by the spectrum of a standard star observed at 
similar air mass. The resulting spectrum is multiplied by a 
solar spectrum to correct for the standard star absorption 
features. The wavelength calibration used OH sky lines 
from the list in Oliva & Origha (1992). 

We extracted spectra for Bl, B2 and C with beam size 
1.2" X 1.2". Since conditions were non-photometric, the 
flux calibration was performed for each source by extract- 
ing photometry with the same beam sizes on the NIC- 
MOS images taken through fihers NIC2 F212N and NIC2 
F222M. The continuum fluxes at 2.142 ^m (for the ob- 
served wavelength of H2) and 2.188 /xm (for the observed 
wavelength of Br7) were interpolated between these two 
filters. The line fluxes were computed by multiplying the 
observed equivalent width of the line by the continuum 
flux. Fluxes and equivalent widths for the Br7 and H2(l- 
O)S(l) emission lines are presented in Table 2. 

Figures 4a and 4b show the fully-reduced spectra. To 
check the NICMOS photometry, we fitted the continuum 
slope between 2.22 /xm and 2.37 /xm using the NIC2 F222M 
and NIC2 F237M fluxes obtained for the corresponding 
apertures. The fits are shown in these figures as dashed 
lines; the agreement is very good. For the ^^C0(2, 0) band 
at 2.293 /xm we measure the spectroscopic index as defined 
in Kleinmann & Hall (1986), which is then transformed to 
a photometric CO index. Measurements of both indices 
for Bl, B2 and C are given in Table 2. Source C does not 
show CO bands. Our values are in good agreement with 
Shier, Rieke, & Rieke (1996) and Vanzi, Alonso-Herrero, 
& Rieke (1998) for the components in common with both 
works. 

2.5. Optical Spectroscopy 

K. Luhman obtained optical spectra of NGC 3690 for 
us with the Red Channel Spectrograph at the Multiple 
Mirror Telescope on Mount Hopkins, on November 27-29 
1997. He used the 270 g mm"! grating (Ablaze = 7300 A) 
to measure spectra from 4600 A to 8600 A. A 2" x 180" slit 
with pixel size 0.6" pixel"-"^ provided a spectral resolution 
of AA = 18 A. To derive the sensitivity fimction of the ar- 
ray he also observed the standard star BD+8 2015. The 
2-D spectra were subtracted from bias-subtracted frames, 
corrected for the sensitivity function and wavelength cali- 
brated using He-Ar-Ne lamp spectra. We did not attempt 
to flux calibrate the data. 

Spectra were extracted centered at B2 with position an- 
gles PA= 120° and PA= 145°, and centered near Bl. Cen- 
tering the slit on Bl with the acquisition camera proved 
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Table 3 
Optical Spectroscopy. 



[Oni]/H/3 


[O i]/Ha 


[Nii]/Hq [Sn]/Ha 


H«/H/3 


EW(H/3) 


EW(Ha) 






Source Bl + B16? Ay 


= 4.2 mag 






2.59 


0.055 


0.20 0.27 


12.4 


15 


170 


Source B2 Ay = 3.5 mag 


1.27 


0.020 


0.34 0.21 


9.8 


5 


55 


Source C Ay = 2.1 mag 


1.42 


0.018 


0.35 0.17 


6.0 


93 


607 


Source C Ay = 2.1 mag 


1.16 


0.042 


0.31 0.29 


6.0 


45 


293 






Region north of A Ay 


= 2.9 mag 






1.76 


0.078 


0.44 0.44 


7.8 


9 


76 






Hll region 6.6" NW of Bl 


Ay = 2.9 mag 






1.57 


0.013 


0.34 0.19 


7.6 


48 


300 






Hll region 7.8" NW of Bl 


Ay =2.1 mag 






1.27 


0.016 


0.47 0.17 


5.9 


79 


550 






Hll region 4.8" NW of B2 


Ay = 1.7 mag 






1.16 


0.016 


0.34 0.21 


5.2 


45 


277 






H II region 2.4" NW of B2 


Ay = 1.9 mag 






0.75 


0.023 


0.31 0.18 


5.5 


36 


243 
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Fig. 4a. — Low resolution i^'-band spectroscopy for Bl, B2 and C. The dashed lines show the Unear fits to the continuum between 2.22 
and 2.37 /im as measured through the NIC2 F222M and NIC2 F237M filters. 




2.3 



2.32 



2.34 
X (/Lim) 



2.36 



Fig. 4b. — Same as Figure 4a but for high resolution JC-band spectroscopy. 
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Fig. 5. — Optical spectra (in arbitrary units) for Bl, B2, C, and C, along with a number of Hll regions. 
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Fig. 5. — Continued. 
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difficult because this nucleus is barely seen in the optical. 
Therefore, the coordinates obtained from the NfCMOS im- 
ages were used, and the slit was kept at the same position 
angles as for B2. Finally the slit was centered at C with 
position angle PA= 80°. The position angles of the slits 
were chosen because PA= 145° includes Bl and B2, and 
PA= 120° includes the Hii regions located northwest of 
Bl. The spectrum centered at Bl probably contains the 
source B16 which is located ~ 1.5" southwest of Bl (see 
Section 3.1). The spectrum of B2 was extracted at the 
position of the peak of the continuum near Ha. The spec- 
trum centered at C at PA = 80° also contains C and 
regions southwest of the nucleus of IC 694. 

For extracting the 1-D spectra we used 3 pixel (1.8") 
apertures. The positions of the slits and the regions for 
which we extracted spectra are shown in Figure 5. The 
fully reduced optical spectra (not calibrated in flux) are 
shown in Figures 5. a and 5.b for those regions for which 
we present line ratios and equivalent widths in Table 3. 
The optical line ratios of all the regions are consistent with 
those of H II as plotted in diagnostic diagrams such those 
in Ho, Filippenko, & Sargent al. (1993). The visual ex- 
tinctions were computed from the Balmer decrement. 

3. RESULTS 

3.1. Continuum Morphology 

Our NICMOS images of Arp 299 reveal an extremely 
complex morphology. Over the past few years it has been 
claimed that the nucleus of IC 694 could harbor an ob- 
scured AGN. An important result of the NICMOS imaging 
is that no dominant point-like source is seen anywhere in 
the system that would be a logical AGN candidate. The 
NICl FllOM image (Figure 2.1) of IC 694 shows that the 
nucleus is resolved into at least three individual sources 
within a region of about 0.86" (~ 175 pc) in diameter. 
These sources could be either young stellar clusters, or 
extinction shadows against a smooth source. Source Bl 
(see Figure 7b for identifications) in NGC 3690 appears 
to be extended at our resolution in the NICl FllOM im- 
age, as docs source B2. The curves of growth for sources 
Bl, B2 and B21 for filter NICl FllOM are shown in Fig- 
ure 6 to illustrate that B2 is clearly extended, whereas 
Bl appears marginally extended. As comparison we plot 
the curve of growth of the source B21 which appears as a 
point source at our resolution. Sources Bl and B2 are also 
surrounded by a number of faint sources (see Figure 7b 
for identifications and Table 4 for the relative distances to 
either Bl or B2). Source C (which is seen as one source 
from the ground) is now resolved into 3 sources that we 
will refer to as CI, C2 and C3 (sec Table 4 for relative 
positions). CI and C2 are extended whereas C3 seems 
to be point-like at the resolution of the infrared images. 
The optical WFPC2 image, because of its higher spatial 
resolution, shows however that C3 is also extended. Ad- 
ditional sources are located southwest of the C complex 
(i.e., C4 and C5) and southeast of Bl (sources Dl, D2, 
D3). The galaxy IC 694 also shows faint sources, most of 
them located along the spiral arms (probably H ii regions), 
and labelled as Al, A2, A3, A4, A5 and A6 (Figure 7a). 
We can place an upper limit to the size of the unresolved 
sources from the NIC FllOM image of 0.08" (~ diffrac- 
tion limit) or ~ 18 pc for the assumed distance. However, 



some of these faint sources appear to be slightly resolved 
in the WFPC2 imagv. 

A comparison between the NICMOS images and the op- 
tical WFPC2 F606W image shows dramatically the effects 
of extinction. The latter image contains the Ha emission 
line, so that low-extinction bright H ii regions will show up 
in this filter. The most remarkable difference appears in 
the component IC 694. A careful alignment of the optical 
and near-infrared images to determine the position of the 
near-infrared peak on the optical image of IC 694 shows 
that the source (presumably the nucleus of the galaxy) is 
completely obscured in the optical (see Figure 2. a). The 
spiral arm southwest of A is clearly detected in the op- 
tical image. The WFPC2 image of NGC 3690 shows a 
closer correspondence with the near-infrared image than 
for IC 694, but still some differences are present. The very 
high spatial resolution of the WFPC2 image shows that 
both sources Bl and B2 are detected and extended (note 
the morphology of the ground-based continuum image near 
Ha where the brightest source is B2, Figure 3), although 
B2 is much brighter than Bl at optical wavelengths (as 
already known from ground-based observations). Most of 
the faint sources seen in the NICMOS images (i.e., the 
B's, C's and D's sources) are also seen in the optical im- 
age, however some of them must be very obscured because 
they are very faint in the optical. On the other hand, 
source B16, which is barely seen in the NIC2 F160W im- 
age, appears to be almost as bright as Bl in the optical. 

The effects of extinction are even more dramatic when 
the NICMOS images are compared with the ultraviolet 
FOG F220W image. We carefully cross-correlated the 
UV sources reported by Meurer et al. (1995, see their 
table 10), with the infrared ones, using relative coordi- 
nates. The brightest source in the UV image (NGC3690-1 
in their notation) is identified with our source B21. The 
UV source NGC 3690-2 corresponds to B2. The double 
source NGC 3690-5 and NGC 3690-8 is located at the po- 
sition of source B16, and therefore does not correspond to 
source Bl (as Meurer et al. 1995 indicated). Source Bl ap- 
pears totally obscured at ultraviolet wavelengths. The po- 
sitions of sources NGC 3690-6, NGC 3690-7, NGC 3690-9 
and NGC 3690-10 are approximately consistent with those 
of C2, C3, C5 and C4 in our notation. Also noticeable 
is that the D complex is relatively bright in the UV, al- 
though it does not stand out in the IR. All the UV sources 
detected in NGC 3690 are embedded in a UV background 
which extends to greater radii than the individual sources. 
IC 694 is very faint in the ultraviolet image (not shown 
here) with only the southeast spiral arm and some very 
diffuse and faint emission near A. 

These comparisons are a warning about the use of the 
UV to study starbursts, particularly ones that are less well 
resolved than Arp 299. The UV light can be dominated by 
lightly obscured regions that represent only a fraction of 
the true starburst region and which may be atypical of the 
integrated properties (e.g., the extended UV backgrounds 
in NGC 3690 and the relative suppression of IC 694). 

3.2. Emission Line Morphology 

The [Fell] 1.644 /im and the Paa continuum-subtracted 
images show bright emission from A, Bl, and C (Fig- 
ures 2. a and 2.b). From both these images it appears that 
very little line emission originates from source B2, confirm- 
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Fig. 6. — Curves of growth for sources Bl, B2 and B21 for the NICl FllOM filter. 
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Fig. 7a. — Contour plot on a linear scale of the NIC2 F160W image of IC 694. We mark the positions of the A sources (Table 4). Orientation 
is north up, east to the left. The field of view is 19.5" x 19.5". 
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Table 4 

Photometry of A, Bl, B2, CI, and faint sources. Relative distances and fluxes. 



Source 


ARA 
s 


ADec 
II 


/(FllOM) 
mjy 


/(F160W) 
mJy 


/(F222M) 
mJy 


/(F237M) 
mJy 


/(F606W)^ 
erg cm~^ A 










0.886 


3.504 


7.280 


7.458 




Al 


0.17 


-1.62 


0.036 


0.100 


0.159 


0.121 




A2 


-0.14 


-2.39 


0.057 


0.124 


0.129 


0.097 




A3 


0.41 


-5.70 


0.047 


0.086 


0.097 


0.075 




A4 


0.48 


-6.20 


0.023 


0.063 


0.073 


0.058 




A5 


-0.27 


-4.71 


0.061 


0.098 


0.108 


0.085 




A6 






0.050 


0.079 


0.080 


0.057 




Bl 


0. 





0.824 


2.909 


10.739 


13.827 


0.19 X 10-15 


Bll 


0.04 


-1.00 


0.158 


0.325 


0.362 


0.333 


0.24 X 10-16 


B12 


-0.08 


-0.42 


0.072 


0.279 


0.372 






B13 


-0.05 


0.91 


0.087 


0.254 


0.397 


0.298 




B14 


-0.11 


0.64 


0.060 


0.104 








B15 


0.09 


0.83 


0.044 


0.115 


0.138 






B16 


-0.06 


-1.10 


0.086 


0.061 






0.10 X 10-15 


B2 








3.678 


5.737 


6.676 


5.525 


0.89 X 10-15 


B21 


0.07 


2.01 


0.440 


0.787 


0.806 


0.663 


0.23 X 10-15 


B22 


-0.17 


-0.85 


0.123 


0.274 


0.338 


0.252 


0.90 X lO-i'' 


B23 


-0.03 


1.14 


0.088 


0.173 


0.152 


0.127 


0.18 X 10-16 


B24 


-0.07 


-0.77 


0.054 


0.091 


0.093 






B25 


-0.30 


-0.60 


0.042 


0.096 


0.104 




0.16 X 10-16 


B26 


0.16 


3.10 


0.039 


0.067 


0.060 




0.20 X 10-16 


CI 








0.199 


1.460 


8.4071 


11.152 


0.72 X 10-16 


C2 


-0.05 


-0.31 


0.171 


0.255 






0.22 X 10-15 


C3 


-0.16 


-1.00 


0.199 


0.307 


0.307 


0.206 


0.76 X 10-16 


C4 


-0.20 


-3.25 


0.119 


0.224 


0.208 


0.154 


0.51 X 10-16 


C5 


-0.23 


-2.63 


0.051 


0.084 


0.083 




0.32 X 10-16 


C 


0.68 


0.77 




0.081 


0.165 


0.160 




Dl 


-1.09 


-0.66 


0.087 


0.136 


0.134 


0.087 


0.47 X 10-16 


D2 


-0.96 


-0.31 


0.013 


0.030 


0.027 






D3 


-0.91 


-0.94 


0.031 


0.048 


0.046 




0.10 X 10-16 



Note. — Relative distances for Al through A6 are from A, relative distances for Bll through B16, and Dl through D3 are from Bl, 
relative distances for B21 through B25 are from B2 and relative distances for C2 through C5 are from CI. The photometry for A, Bl and 
B2 is through a 0.9"-diameter photometry, whereas that of the faint sources is for a 0.5"-diameter aperture (corrected for aperture effects). 
Errors for the measured fluxes are 7% for the bright sources and 10% for the faint sources. 
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Table 5 





Observed iJ-BAND magnitudes, 


COLORS AND 


(CO)nicmos index. 




Source 


H 


J -H 


H-K 




CO 


R-J 


A 


13.724 


1.330 


1.361 




-0.025 




Al 


I 7 ^^iS 

I I .ooo 


±.U 1 o 


1.039 




u.zoy 




A2 


17.357 


0.908 


0.570 




0.249 




A3 


17.746 


0.807 


0.649 




0.229 




A4 


18.095 


1.049 


0.686 




0.203 




A5 


17.614 


0.720 


0.632 




0.209 




A6 


17.846 


0.701 


0.539 




0.256 




Bl 


13.927 


1.230 


2.061 




-0.288 


2.18 


Bll 


16.305 


0.880 


0.638 




0.079 


2.38 


B12 


16.473 


1.291 


0.837 








B13 


16.575 


1.108 


1.019 




0.249 




B14 


17.539 


0.770 










B15 


17.430 


1.034 


0.720 








B16 


18.121 


0.190 








-0.15 


B2 


13.190 


0.699 


0.687 




0.172 


1.84 


B21 


15.347 


0.788 


0.548 




0.178 


1.15 


B22 


16.493 


0.930 


0.751 




0.252 


2.97 


B23 


16.991 


0.853 


0.386 




0.165 


2.09 


B24 


17.687 


0.750 


0.541 








B25 


17.636 


0.939 


0.618 






1.57 


B26 


18.028 


0.766 


0.418 






1.08 


CI 


14.676 


1.705 


2.646 




-0.327 


1.97 


C2 


16.570 


0.672 








0.13 


C3 


16.368 


0.694 


0.522 




0.331 


1.37 


C4 


16.712 


0.823 


0.443 




0.259 


1.34 


C5 


17.781 


0.727 


0.510 






0.89 


Dl 


17.253 


0.703 


0.505 




0.352 


1.04 


D2 


18.877 


0.949 


0.392 








D3 


18.372 


0.699 


0.462 






1.55 


Note.— 


-Errors in the colors 


are ±0.10 and ±0.14 for the bright and faint 


sources respectively, whereas 



the errors associated with the (CO)nicmos are of the order of 20%. 
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Fig. 7b. — Contour plot on a linear scale of the NIC2 F160W image of NGC 3690. We mark the positions of the B, C and D sources 
(Table 4). Orientation and field of view as in Figure 7a. 



ing the result of Fischer, Smith & Glaccum (1991). The 
[Fe ii] 1.644 /.mi and Paa images trace the emission from 
H II regions located in the spiral arms SE of IC 694. The 
Paa emission of the nucleus of IC 694 is similar to the ra- 
dio H92a morphology (Zhao et al. 1997). However, the SE 
emission seen in Zhao et al. (1997) map does not correlate 
exactly with the Paa morphology, possibly due to some 
obscuration effects at the infrared wavelengths. There are 
a number of Hii regions northwest of Bl, and east of C 
reaching all the way to source C and extending toward 
IC 694. In the Paa image we can also resolve the B16 
source. 

The H2 image of IC 694 shows intense emission from the 
nucleus of the galaxy together with a beautiful butterfly- 
like diffuse emission extending over approximately 3 arc- 
sec (see Figure 2. a). The molecular hydrogen image of 
NGC 3690 shows mainly point-like emission from sources 
Bl, B2 and C, and some diffuse emission from source C 
(Figure 2.b). 

3.3. Photometry of A, Bl, B2, C, and the faint sources 

To study the faint sources in Arp 299, we need accurate 
magnitudes and colors. However, some of these sources, 
such as the faint objects near Bl and B2, are located on 
a local background with a substantial gradient. To sub- 
tract the contribution of the underlying galaxy, first we ex- 
tracted subimages of each source, then we produced a two 
dimensional model of the local background for each source. 
The source plus background fluxes were fitted unidimen- 
sionally along both rows and columns. The source was rep- 
resented with two Gaussians with different FWHMs and 
amplitudes, whereas we used a linear fit to the local back- 

V2 = (■^'="duai(ij)-<residual»^ n = i X j and residual = i: 



ground. Although the NICMOS PSF is not perfectly rep- 
resented by Gaussian functions, the one-dimensional emis- 
sion profiles were found to be fitted well. The two fits along 
rows and columns were then averaged together to produce 
independent 2-D source and the background models. The 
model background was subsequently subtracted from the 
original sub-image. 

Aperture photometry was done on the background sub- 
tracted images. Most of the faint sources appear to be 
unresolved, so the photometry was obtained through rel- 
atively small apertures (typically 0.5" in diameter). The 
final magnitudes and fiuxes were corrected for the behav- 
ior of the PSF so that the measurement included all the 
flux associated with a point source. The errors in the pho- 
tometry for each individual source were computed as the 
standard deviation of the subimage resulting from the sub- 
traction of the background+source model image from the 
original image^. These errors are approximately 0.1 mag 
for the faintest sources, decreasing to 0.07 mag for bright 
ones. 

Fluxes (in mJy) measured through the NICl FllOM, 
NIC2 F160W, NIC2 F222M and NIC2 F237M fifiers along 
with the fluxes (in erg cm~^ s~^) through the WFPC2 
F606W filter are given in Table 4. For the WFPC2 F606W 
filter we used the photometric calibration from the header 
of the image, since this filter was not included in the set 
of WFPC2 cahbrated filters. In Table 5 we give the H- 
band magnitudes and J — H, H — K and R — J colors. 
The J — H and H — K colors were obtained by convolving 
an empirical galaxy spectral energy distribution with both 
the NICMOS and J, H, K filter transmission functions and 
detector response functions. The same convolutions were 

e - model 
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Table 6 

Extinction to the gas for the different components of Arp 299. 



Aperture 


source 


/(Paa) 


/(Paa)//(Ha) 


Ay (screen) 






erg cm~-^ s~^ 




mag 


circular 5" 


A 


3.7 X 10-13 


1.9 


5 




Bl 


2.3 X 10-13 


0.5 


2 




C 


4.4 X 10-13 


0.4 


2 


4.5" X 2" 


A 


3.1 X 10-1^ 


1.5 - 3.0 


4-5 






1 7 V 1 fl— 13 
1. 1 X lU 


U.o — l.U 


O — 4 






3.3 X 10 


0.4 — 0.8 


2 — 3 


2" X 2" 


A 


2.4 X 10-1^ 


3-5 


5-6 




Bl 


1.3 X 10-13 


0.7-0.9 


3-4 




C 


2.7 X 10-13 


0.6 - 0.9 


3-4 


Aperture 


source 


/(Paa) 


/(Br7)//(Paa) 


Ay (screen) 






erg cm ^ s l 




mag 


2.4" X 4.7"'' 


A 


3.5 X 10-^3 


0.115 


24 ±7 


1.2" X 1.2" 


Bl 


8.3 X 10-14 


0.102 


16 ±7 




C 


1.8 X 10-13 


0.142 


37 ± 10 


1.1" X 2.2"a 


Bl 


9.6 X lO-i"! 


0.104 


17 ± 7 


Aperture 


source 


/(Paa) 


/(Paa)//(Pa/^) 


Ay (screen) 






erg cm-^ s-l 




mag 


2.4" X 4.7"" 


A 


3.5 X 10-13 


5.7 


9±2 



'^Fltixes of Br7 calculated from EW of Br7 and NICMOS continuum calibra- 
tion (Craig Kulesa private communication). 

''Flux of Pa/3 calculated from EW of Pa/3 and NICMOS continuum calibra- 
tion (Craig Kulesa private communication). 



also computed for the galaxy reddened by as much as 40 
magnitudes. The convolutions were compared to the val- 
ues for the standard star P330E in each filter set to derive 
a transformation between the NICMOS F110M-F160W- 
F222M system and the JHK system. This procedure also 
corrected for the large color difference between the stan- 
dard star and our objects in the rather broad F160W filter. 
The Johnson R magnitudes are interpolated between the 
fluxes at 6060 A and at 1.10 /im. 

In addition we measure a photometric NICMOS CO in- 
dex defined as, 

/(F222M) - /(F237M) 
COnicmos = ;(F222M) 

This quantity does not formally correspond with the stan- 
dard photometric CO index, but as Figure 8 illustrates, it 
is nearly equal to the photometric CO index. The relation 
plotted in Figure 8 was derived by integrating the SED 
for an old stellar population across both the two NICMOS 
filter functions and the ground-based filter functions. The 
SED was redshifted to .010 and also reddened by amounts 
ranging from Ay = mag to Ay = 40 mag. The measured 
values for both the bright and faint sources are given in 
Table 5. The NICMOS CO index is subject to the same 
caveats as the ground-based index because it measures the 
slope of the continuum between 2.22 /tm and 2.37 /xm as 
well as the strength of the CO absorption. As we will 
discuss in Section 3.3.2 it is sensitive to both extinction 
effects and dilution of the CO bands by hot dust emission. 



3.4. Extinction 

To interpret the observations requires an accurate es- 
timate of the extinction to each source. The extinction 
to the gas can be estimated from hydrogen recombination 
line ratios, so long as the comparison is based on measure- 
ments with the same beam sizes. To determine the ex- 
tinction to the stars one can compare observed colors with 
the colors of an evolved stellar population. The infrared 
colors of evolved stellar populations are nearly universal, 
independent of age. In addition a correction for the dust 
contribution to the /T-band fluxes -if any- is required prior 
to dcreddening the stellar colors. 

The largest uncertainty in the extinction by far comes 
from the unknown distribution of dust within the source. 
Meurer et al. (1995) from their study of the UV proper- 
ties of starburst galaxies concluded that the dust geome- 
try has to be in a foreground screen configuration near the 
starburst. However, UV-selected galaxies will be biased 
toward low extinction and possibly an unrepresentative 
distribution of the dust. In addition, as emphasized in Sec- 
tion 3.1, the UV may be dominated by non-representative 
regions within a galaxy. Furthermore, within this study 
NGC 3690 was a notable exception because it appeared to 
have patchy extinction. Other authors (e.g., McLeod et 
al. 1993 for M82; Genzel et al. 1995 for NGC 7469; and 
Sugai et al. 1999 for Arp 299) find that a simple foreground 
screen model does not provide a good fit for relatively dust 
embedded starbursts. They used a model in which dust 
and gas are homogeneously mixed, and showed that the 
foreground screen model gives a lower limit - sometimes 
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Fig. 8. — Relation between the photometric (CO) index and the (CO)nicmos for an elliptical galaxy redshifted to z = 0.01 and different 
values of the extinction from Ay = mag through Ay = 40 mag. 



Table 7 

Extinction to the stars, corrected colors and NICMOS CO 
INDICES. Dust emission. 



Source 


Av 


(11 - A-)c„rr 


CO^^I( 'M( )s (corrected) 


(/<luj,t /,/t.o(al)/\' 


(1) 




(■'>) 


(1) 


(•') 


A 


6.5 


1.01 


0.03 




Bl 


5.2 


1.78 


-0.24 




CI 


9.9 


1.90 


-0.23 




A 


24 ±7 


0.08 ± 0.40 


0.15 ± 0.04 





Bl 


17 ±7 


1.20 ± 0.40 


-0.15 ±0.04 


~ 0.6 


B2 


3.5 ±0.5 


0.31 ± 0.03 


0.19 ±0.01 





CI 


37± 10 


0.43 ± 0.70 


0.00 ± 0.06 


< 0.2 




Note. — The errors in columns (3) and (4) take into account the uncertainties as- 
sociated with the extinction. 
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by a large factor - to the extinction. 

3.4.1. Extinction to the gas 

We use a foreground dust screen model and consider its 
effects on the hydrogen recombination lines. For the in- 
frared extinction law we use a polynomial fit to the Aj, 
Ah and Ak values given in He et al. (1995). The re- 
sulting relative extinctions are, A(Pa/3) = 0.254 x Ay, 
A{Paa) = 0.128 X Ay, AiBi-/) = 0.114 x Ay. For the opti- 
cal lines we use A(H/3) = 1.20xAy and A(Ha) = 0.73xAy 
(Riekc & Lebofsky 1985). The intrinsic values for the Case 
B hydrogen line ratios are taken from Hummer & Storey 
(1987) for physical conditions Tg = 10, 000 K and = 100 
cm~^. 

The results for different lino ratios and aperture sizes 
are presented in Table 6. In this table the fluxes of the 
hydrogen recombination lines are ratioed to the Paa flux, 
which was measured through the corresponding aperture 
from our NICMOS images. The two values of the Paa to 
Ha ratios and extinctions for the 4.5" x 2" and 2" x 2" 
apertures correspond to the uncertainty of the Ha flux cal- 
ibration, as discussed in Section 2.3. The results from the 
Br7 to Paa line ratios arc presented in the second part 
of this table. We did not obtain K-hand spectroscopy for 
source A, but we made use of the data in Kulesa et al. 
(1999). Instead of using their flux calibration, we mea- 
sured the continuum fluxes from the NICMOS images and 
used their equivalent widths of Br7 and Pa/3, to compute 
the line fluxes. The errors listed in Table 6 for the Ay 
derived from the Br7 to Paa and Paa to Pa/3 line ratios 
are simply the propagation of the errors in the measured 
equivalent widths. In Table 3 we give the extinctions de- 
rived from the Balmer decrement. 

From Table 6 the extinction tends to be higher when 
line ratios involving longer wavelengths are employed. 
For instance, the extinction towards source A increases 
from Ay = 5.8 mag from the Balmer decrement up to 
Ay = 24 ± 7 mag for ratios involving longer wavelengths. 
A similar trend is found for source C (Cl-f-C2) where ex- 
tinctions range from Ay = 2.1 mag from the optical line 
ratios (sec Table 3) to Ay = 37 ± 10 mag from the Br7 to 
Paa line ratio. The large error of the latter value is due 
to the small value of the differential extinction between 
the two lines together with a 10% uncertainty in the flux 
calibration. Prom the fluxes of Brackett-series lines given 
in Sugai et al. (1999) we obtain values of the extinction to 
source C of between Ay = 8 mag and Ay = 12 mag. For 
source Bl the extinction does not seem to vary in such 
a dramatic way, but it still increases toward longer wave- 
lengths, with a maximum value of Ay = 17 mag. 

The differing values as a function of wavelength are the 
usual indication that the extinction is either patchy or 
mixed with the emitting sources along the line of sight (or 
both). In confirmation, for C1+C2 patchy extinction is 
seen from the optical/near infrared color map (see next 
section). In addition, our data indicate that the extinc- 
tion increases for many sources as the aperture is reduced 
in size. In such a highly-obscured and complex environ- 
ment we are not probing the same regions in the optical 
and infrared even when the same aperture sizes are used. 
This conclusion was already implicit in the differing mor- 
phology of the continuum images between these spectral 
regions (Section 3.1). 



The relative distribution of gas and dust appears to be 
similar for both molecular and ionized hydrogen. Our es- 
timates from near infrared hydrogen recombination lines 
and a foreground screen are in good agreement with the 
values derived from the near infrared molecular hydrogen 
lines and a foreground screen by Kulesa et al. (1999) which 
are: Ay = 23 mag for source A and Ay = 15 mag for both 
sources Bl and C. 

3.4.2. Extinction to the stars 

We can obtain a qualitative idea about the differen- 
tial extinction from the color maps presented in Figure 9 
(infrared-infrared H — K and optical-infrared Mpeoew — 
H). Because the WFPC2 F606W filter contains Ha, the 
optical to near-infrared color maps will show in great con- 
trast bright unobscured H ii regions. Both color maps con- 
firm that the extinction to the nucleus of IC 694 is very 
high. It is interesting that the H — K map shows a similar 
morphology to the molecular hydrogen emission (the but- 
terfly pattern). The Mpeoew ^ H color map of region C 
shows patchy extinction with a very red embedded source 
surrounded by bright Hii regions. The ground-based Ha 
image shows very bright emission from region C, which 
suggests that the Ha emission is dominated by the bright 
outer Hii regions, whereas the near-infrared recombina- 
tion lines trace the emission from the very red embedded 
source. Near Bl we can clearly see the bright Hll re- 
gion associated with source B16, whereas the other sources 
(Bl, and Bll through B15) seem to be obscured by a dust 
screen, rather than patchy extinction. 

A more quantitative estimate of the extinction to the 
stars can be obtained from the colors. The optical-infrared 
color R — J is sensitive to age affects and the sources 
in Arp 299 appear to be young. We therefore use the 
J — H and H — K colors, which are relatively unaffected 
by age. As found for the extinction to the gas, when a 
simple foreground screen model is assumed, the values of 
the extinction increase from the J — H values to the longer 
wavelength H — K ones. This effect is shown in the first 
part of Table 7, where we give the extinction inferred from 
dereddening J — H to the color of a normal stellar popu- 
lation {J — H = 0.7) along with the corrected H — K and 
(CO)nicmos values using that extinction. For the F237M 
fluxes we interpolated the extinction between the values 
for Ak and A^ given in Rieke & Lebofsky (1985). The 
dereddened H — K colors are still very red, but the ques- 
tion remains whether this is the effect of hot dust emission 
and/or differential extinction. In the second part of Ta- 
ble 7 we give the same quantities but corrected for the 
highest possible extinction derived from the hydrogen re- 
combination lines. 

We can distinguish between these extreme values of the 
extinction by computing (CO)nicmos versus H — K as 
a function of the properties of any emission by hot dust. 
Figure 10 shows the results for sources A, Bl, B2 and CI. 
The open star symbols are the values obtained for the ex- 
tinction derived from the observed J — H color, whereas 
the fllled star symbols are the values obtained using the 
highest extinction derived from the hydrogen recombina- 
tion lines. The solid lines represent mixing curves of a 
normal stellar population and hot dust with different tem- 
peratures. For the normal stellar population we integrate 
the spectrum of an elliptical galaxy redshifted to ^ = 0.01 
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using the filter transmission curves and detector response 
functions. 

The H — K and (CO)nicmos values were obtained com- 
paring to the values for the standard star P330E as de- 
scribed in Section 3.3. For the redshifted elliptical galaxy 
wc obtain H - K ^ 0.23 and (CO)nicmos = 0.13. The 
dashed lines are the locii of the ratio (/dust//totai)if with 
values 0.2, 0.4, 0.6, 0.8 and 1. These mixing curves are 
computed by adding different dust contributions (in the 
form of a blackbody) to the elliptical galaxy. Figure 10 
shows that reproducing the (CO)nicmos indices and the 
H — K colors with the extinction obtained by dereddening 
the J — H colors would require dust temperatures higher 
than T = 1000 K for all sources. Such behavior is atypical 
of starbursts (e.g., Rieke et al. 1993). 

An independent approach to the same issue can be found 
in Satyapal et al. (1999). They use measurements of the 
3.29/im PAH emission feature and a standard ratio of fea- 
ture strength to hot dust continuum to conclude that hot 
dust can contribute only 7, 3, and 12% respectively to the 
K-band fluxes of sources A, B, and C. 

The other possibility is that there is differential extinc- 
tion to the stars, as we have already found for the gas. 
Using the highest values for the gas extinctions (Table 6), 
we obtain more reasonable temperatures for the dust to- 
gether with a smaller dust contribution to the i^-band 
flux. Indeed, only source Bl needs a dust contribution 
of the order of 40%, whereas the other sources require 
(/dust//totai)K < 0.2. The fact that using the highest ex- 
tinctions derived for the gas does produce reasonable val- 
ues for the dereddened values oi H — K and (CO)nicmos 
suggests that the extinction to the gas and to the stars 
is similar (at least in the near- infrared) . The general ten- 
dency for the extinction in the gas to be about two times 
higher than that to the stars (e.g., Calzetti, Kinney, & 
Storchi-Bergmann 1994) does not appear to hold for this 
galaxy. 

A more realistic approach to the treatment of the ex- 
tinction should include a model in which stars and dust 
are mixed (e.g., Witt, Thronson, & Capuano 1992; Witt 
& Gordon 1999). We used the recently published mod- 
els of Witt & Gordon (1999) to evaluate the extinction 
to source A. We find that the H — K color excess can be 
only reproduced with their shell models, using a K-h&nd 
extinction of Ak — 1.8. This is equivalent to the value 
used to dercdden the /C-band magnitude of component A 
when using the dust screen model (see Section 4.5), con- 
firming the assumption that the system becomes roughly 
optically thin near the iiT-band. Note that in Witt & Gor- 
don (1999) models, the relation between Ak and Ay is 
model dependent. 

Finally since (/dust/ZtotaO^: for source C appears to be 
relatively small, the lack of CO bands is not due to dilu- 
tion by hot dust emission, but rather due to the youth of 
the stellar population that dominates the /T-band light. 
This conclusion is confirmed by the lack of CO bands in 
the if-band also (see spectrum in Vanzi et al. 1998). 

4. STAR FORMATION PROPERTIES 

We can distinguish at least six star forming environ- 
ments in the Arp 299 system: 1.) Component A, the nu- 

^<f>{m)Am oc m~^'^^dm between 1 and 80 Mq, and 0(m)dm oc m~ 



cleus of IC 694; 2.) component Bl; 3.) B2; 4.) C; 5.) C; 
and 6.) additional bright Hii regions identified individu- 
ally in our Paa images. In this section we will make use 
of evolutionary synthesis models (Rieke et al. 1993; En- 
gelbracht et al. 1996, 1998) and other arguments to study 
the star formation properties of these regions. We will as- 
sume a truncated Salpeter IMF"^. Normalized to the same 
total mass, this IMF is virtually identical to IMF8 found 
by Rieke et al. (1993) to give a good fit to the starburst 
properties of M82. We first consider the parameters that 
allow us to deduce the star forming properties of these 
sources. 

4.1. Luminosity 

Since Arp 299 is an infrared luminous galaxy, we can 

assume that most of the bolometric luminosity is emit- 
ted in the infrared. An estimate of the fraction of the 
total infrared luminosity from each component can be ob- 
tained from ground-based mid-infrared observations of the 
system (Gehrz et al. 1983; Wynn- Williams et al. 1991; 
Keto et al. 1996). The relative fluxes at 10/im, 12 /im, 
20 /im and 32 /im through a 6"-diameter aperture for 
each component arc: A = 27% 35% 42% ^ 63%, 
Bl = 55% ^ 53% ^ 40% ^ 30%, and C + C = 17% 
12% 18% 7%. We show below that about 10% of the 
total ionizing luminosity is generated in supergiant H II 
regions that would not be detectable individually in the 
ground-based data. Hence we can assume that ~ 50% of 
the total infrared luminosity of the system is concentrated 
at A (a similar fraction was derived by Joy et al. 1989 for 
mid- infrared and IRAS far-infrared data), ~ 27% is con- 
centrated at component Bl {+ B2) of NGC 3690, ~ 13% 
is at C-l-C' area, and ~ 10% in the supergiant Hit regions. 
Prom the mid-infrared maps of Keto et al. (1996), it can 
be inferred that most of the emission for the area B1+B2 
comes actually from component Bl. For our model fltting 
we will take the bolometric luminosities given in Table 8 
as lower limits, to allow for any escaping energy. 

4.2. Mass 

An estimate of the available mass is useful as an input 
to the synthesis models. Shier et al. (1996) derived a 
mass of 5.6 ± 1.8 x IO^Mq for IC 694 from the velocity 
dispersion measured for the 2.3/im CO band. This mass is 
in reasonable agreement (but is larger than) the determi- 
nation from radio emission lines (Zhao et al. 1997). For 
component B2 of NGC 3690 the dynamical mass inferred 
by Shier et al. (1996) is 0.6 x 10^ M©, whereas the shallow 
CO index of Bl did not allow them to get a reliable esti- 
mate for this component. We face a similar problem with 
source C, with its lack of CO bands. 

Starburst luminosity varies more slowly than other pa- 
rameters such as ionizing flux. Since it appears that some 
of the star forming regions in this system are not in nu- 
clear potential wells, the only useful mass associated with 
them is the actual mass associated with the local star for- 
mation. To estimate the masses of Bl and C where there 
are no dynamical measurements, we assume the mass-to- 
light ratio derived for the starburst model of component A, 
M/LiR = 0.004. The masses of each component are given 
in Table 8. For Bl we subtracted the dynamical mass of 

^dm between 0.1 and 1 Mq. 
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H-K 



Fig. 10. — (CO)nicmos versus H — K . The solid lines represent mixing curves of hot dust and a normal stellar population rcdsliiftcd 
to 2 = 0.010 (H — K = 0.23 and (CO)nicmos = 0.13), whereas the dashed lines are the fraction of hot dust emission in the K-hsa\d: 
{ fduat / f total) K = 0.2, 0.4, 0.6, 0.8 and 1 (from bottom to top). The open stars are the values for H — K and (CO)nicmos corrected for 
the extinction derived from the J — H color, whereas the filled stars are the values tor H — K and (CO)nicmos corrected for the (highest) 
extinction derived from the Br7/Paa line ratio. 



B2. 

The mass of molecular gas in the interacting system 
Arp 299 has been estimated from radio CO observations 
assuming a standard Galactic CO to H2 conversion (see 
for instance Sargent & Scoville 1991). These estimates 
obtain the troubling result that the mass of molecular 
gas accounts for most of the dynamical mass (Mgas — 
0.7 X Mdyn). As discussed by Shier et al. (1994) the stan- 
dard conversion factor is probably a factor of 3 to 10 too 
high for infrared luminous galaxies. The work of Shier 
ct al. (1994) has been questioned by Scoville, Ytm, & 
Bryant (1997) and Solomon et al. (1997) on the basis that 
2 (Um spectroscopy will not penetrate the interstellar dust 
in these systems and hence will underestimate the dynam- 
ical mass. However, Zhao et al. (1997) used the H92a 
lino to meastiro a dynamical mass for IC 694 loss than, 
but consistent with, that obtained by Shier et al. (1994), 
supporting the general validity of the arguments of Shier 
et al. (1994). Braine & Dumke (1998) have recently esti- 
mated H2 masses from the mm-wave thermal emission of 
interstellar dust in Arp 299 and similar systems and find, 
in agreement with Shier et al., that the true H2 masses 
are probably factors of 3 to 10 lower than given by the 
standard CO conversion. We therefore assume for all the 
nuclei that the conversion to M(H2) is 30% of the standard 
value. The molecular mass in A and probably B2 is then 
about 15% of the total mass, consistent with theoretical 
expectations for when molecular gas in a galaxy nucleus 
would tend to break up into star forming clouds (Wada & 
Habe 1992, Bekki 1995). 

4.3. The Number of Ionizing Photons 



Zhao ot al. (1997) used observations of the radio hy- 
drogen recombination lino H92q; to model the physical 
conditions of the H II regions in components A and Bl of 
Arp 299. Their modeling has a number of free parameters 
(such as size of the Hll region, electron density, etc). The 
advantage in using this radio line is that it is not affected 
by extinction, although the estimated number of ionizing 
photons is slightly dependent on the model parameters 
(see Zhao et al. 1997 for a discussion). The resulting un- 
certainties do not allow them to distinguish between two 
oxtromo models involving either a small number (a few 
tens) of very luminous Hii regions with typical sizes of 
25 pc, or a large number (~ 10^) of modest luminosity H II 
regions of very compact sizes (~ 0.08 pc). 

The values derived by these authors are TVly = 3.8 — 
4.7 X 10^** s^^ for IC 694 (component A). These values cor- 
respond to those from the infrared lines if 16 mag < Ay < 
18 mag. For component Bl of NGC 3690, the fits of Zhao 
ct al. yield Ni^y = 1.5 - 8.8 x lO-'^'^s-^ The larger value 
holds for a small number of large H 11 regions, which is as 
indicated by our Paa imaging and is particularly likely for 
source Bl because of its compact size in our images. It 
yields an estimate of Ay ~ 11 mag. 

Within the errors, these extinction values agree with 
those obtained from the longer wavelength infrared lines. 
They are significantly larger than would have been de- 
duced from the Balmer decrement or even from the shorter 
wavelength infrared lines such as Pa/3 (A = 1.28 /im). We 
adopt values consistent with both the radio and infrared 
determinations: Ay = 17 mag for A and Ay = 11 mag for 
Bl. Unfortunately Zhao et al. (1997) did not model com- 
ponent C. From our data and those of Sugai et al. (1999), 
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we will adopt Ay = 15 mag. In general, our arguments 
regarding the extinction behavior have close quantitative 
agreement with the work of Sugai et al. (1999). In Ta- 
ble 8 we give the extinction-corrected values for number 
of ionizing photons (A''Ly). 

An independent test can be made from the free-free ra- 
dio continuum. In the case of M82, Pnxley et al. (1989) 
have compared this continuum with the strength of the 
H53a line, which is much less subject to modeling uncer- 
tainties in converting to TVLy than is H92q:. We therefore 
have confidence that the strength of the free-free contin- 
uum is correctly related to N^j, in M82. By correcting for 
distance, we predict strengths of 14, 3, and 11 mJy for 
the free- free emission at 3mm in components A, B, and 
C, respectively, of Arp 299. The measured fluxes at this 
wavelength are 17 ± 2, 5 ± 2, and 9 ± 2 mJy respectively 
(Aalto et al. 1997). The agreement is excellent. 

For component A, our data, those of Sugai et al. (1999), 
the comparison with Zhao et al. (1997), and the free-free 
radio flux are in excellent agreement but strongly disagree 
with the recent conclusion by Satyapal et al. (1999). The 
latter reference uses an extinction of only Ay = 6 and 
deduces that the Lyman continuum flux is three times 
smaller than obtained by the other four approaches. Satya- 
pal et al. acknowledge the discrepancy with the H92a 
flux, but they adopt the smaller ionizing fluxes suggested 
by their extinction levels for their modeling. Since our 
improved infrared data agree with the higher extinction 
levels, we believe their values are underestimated due to 
their use of Pa/3 in the extinction estimate. Evidently, op- 
tical depth effects at 1.28/im are significant. Satyapal et 
al. adopted Ay = 2 mag for component C, whereas both 
our data and those of Sugai et al. indicate a much larger 
value, Ay > 10 mag. The resulting underestimation of the 
ionizing fluxes undermines the modeling of Satyapal et al. 
(1999) for these components. 

Although for simplicity we (and Satyapal et al.) have 
used a foregroimd screen model for the extinction, the 
agreement with Zhao et al. (1997) is only achieved with 
the 'high' extinction values from the long wavelength re- 
combination lines. This behavior demonstrates that the 
dust and gas (and probably the stars) are mixed, and that 
the system becomes roughly optically thin near the K- 
band. 

4.4. The Supernova Rate 

The supernova rate in starburst galaxies can be esti- 
mated from radio observations (assuming that the non- 
thermal radio flux is synchrotron emission produced 
in shocks associated with supernovae), and from near- 
infrared [Fe ii] emission lines using the calibration obtained 
for the starburst galaxy M82 (Vanzi & Rieke 1997; Alonso- 
Herrero et al. 1999). Using the 20 cm radio data of Gehrz 
et al. (1983) for the apertures which most closely match 
ours, and the relation between 5 GHz flux density and su- 
pernova rate measured for M82 by Huang et al. (1994), 
we obtain rates of 0.42, 0.11, and 0.06 yr~^ for sources A, 
B, and C respectively. The Condon & Yin (1990) calibra- 
tion for the supernova rate from Galactic supernovae yields 
similar values. However, the radio spectrum of source A 
is too flat near 20 cm to match typical supernovae; we 
believe it is affected by free-free absorption as previously 
found for M82 (Rieke et al. 1980). We use the relative nu- 



clear flux densities at 8.3 GHz from Zhao et al. (1997) to 
correct for this effect. We arrive at a best estimate; of the 
supernova rates of 0.65, 0.11, and 0.06 yr~^, respectively. 
All of these values may be slightly high because free-free 
emission is expected to contribute at the ^ 10-20% level 
to the observed radio fluxes. 

We can also use the newly determined calibration be- 
tween the [Fell] 1.644 /xm luminosity and the supernova 
rate for M82 (Alonso-Herrero et al. 1999), which gives 
about a two times lower supernova rate than the previ- 
ous best calibration (Vanzi & Rieke 1997). Correcting the 
observed [Fell] fluxes for the appropriate extinctions, we 
derive supernova rates of 0.52, 0.06 and 0.06 yr~^ for A, Bl 
and C (extinction ^ = 15 mag) in apertures of 2.4" x 4.6", 
1.2" X 1.2", and 1.2" x 1.2". The agreement with the radio- 
based estimate is quite good; for the model fitting we use 
the radio results, but there would be no change if we had 
used [Fell]. 

4.5. Results from the evolutionary synthesis models 

We can interpret the starburst properties with evolu- 
tionary synthesis models. The observables that we will try 
to reproduce with these models are the bolometric lumi- 
nosity Lboi, the ii'-band luminosity, the number of ionizing 
photons, the CO index, the supernova rate and the EW 
of Paa. The apertures for all the parameters are matched 
to the infrared spectroscopy. As in Rieke et al. (1993) we 
normalize the output of the evolutionary synthesis models 
to the observed values for each component; we will plot the 
evolution in time of the six quantities so that our target 
quantity will be unity. The mass-to-light ratios for both 
i^-band luminosity and the bolometric; luminosity given 
in Engelbracht et al. (1996) are used to estimate the con- 
tributions of the old underlying population to the total 
K band and bolometric luminosities. For simplicity, we 
assume a single burst of star formation with a Gaussian 
FWHM of 5 X 10^ yr. The peak of star formation occurs 
at 5 x 10® yr after the beginning of the burst. 

Source A. 

We assume an extinction of Ay = 17 mag. The K band 
luminosity and the A^Ly are corrected for extinction us- 
ing a foreground dust screen. Because we find that the 
source becomes roughly optically thin at K, this proce- 
dure should give reasonably accurate results. For the CO 
index, we take the value from Shier et al. (1994), which 
is measured spectroscopically and can be converted accu- 
rately to the photometric index used in the models. The 
fit for the observables of source A (Figure 11. a) uses a total 
mass in newly formed stars of 7 x 10* M0 (~ 12% of the 
dynamical mass) and yields an age for the star formation 
episode of 6 — 8 x 10® yr (after the peak of star formation). 
There is a modest discrepancy in the supernova rate. Mea- 
sured from the rising half power of the assumed Gaussian 
in the star formation rate, the age of the burst is then ~ 
11 Myr. 
Source Bl 

For the assumed level of the extinction {Ay = 11 mag), 
the contribution of hot dust emission to the observed K 
luminosity is (/dust //total) if — 0.5, which needs to be sub- 
tracted to get the stellar X-band luminosity. The observed 
CO index and equivalent width of Paa also need to be cor- 
rected for the dust emission contribution to the continuum. 
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Table 8 

Properties of the components of Arp 299. 



Source 


Mass 


Av 


Lm 


ATLy (total) 


K 


CO 


SNr 


EW(Paa) 




(M©) 


(mag) 


(s-i) 


(s-i) 


(mag) 




(yr-i) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A 


5.8 X 10^ 


17 


2.5 X 10" 


4 X 10^* 


-23.87 


0.17 


0.65 


220 


Bl 


> 0.4 X 10^ 


11 


< 1.4 X 10" 


9 X 10^3 


-22.64 


-0.15 


0.12 


100 


B2 


0.6 X 10^ 


3.5 






-21.76 


0.19 






C 


> 0.2 X 10^ 


15 


0.5 X 10" 


3 X 10^4 


-22.71 


0.00 


0.06 


500 


C 


> 0.1 X 10^ 


10 


0.25 X 10" 


3 X 1053 


-18.98 




0.03 


1500 


SG Hii 


~ 0.02 X 10^ 




0.2 X 10" 


> 2.2 X 10^3 











Note. — (1) SG H li represents the sum of the contributions from the 19 most luminous HII regions in the system. No 
extinction corrections have been applied for them. (2) Dynamical mass for A and Bl from Shier et al. (1996). Masses 
for remaining components represent forming stellar mass to produce observed infrared luminosity. (3) Assumed foreground 
extinction level. (4) Estimated infrared (8 — 1000 /Ltm) luminosity for each component. (5) Total number of ionizing photons 
corrected for extinction. (6) Extinction corrected K-hand absolute magnitude for the small apertures (2.4" x 4.6" for Source 
A, 1.2" X 1.2" for Bl and C, 2" x 2" for C'). (7) NICMOS CO index. (8) Supernova rate for small apertures from radio. 
(9) Equivalent width of Paxx for the small apertures. 




age 



(millions) 



Fig. 11a. — Evolutionary synthesis model fit for source A and a truncated Salpeter IMF. The extinction to the gas is assumed to be the same 
as the extinction to the stars, Ay = 17 mag. The curves are as follows. Solid line (declining) the EW of Paa. Solid line Jf-band luminosity. 
Dotted line the bolometric luminosity. The dashed line the supernova rate. The dot-dash line the CO index. The dot-dot-dot-dash line the 
number of ionizing photons. 



Table 9 
Synthesis Fit to C. 





Lm 


A^Ly (total) 


K 


SNr 




(L0) 


(s-l) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


Observations 


0.25 x 10" 


3 X 10^3 


-18.98 


0.03 


Model 


0.3 X 10" 


4 X 10^3 


-19.0 


0.025 
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Fig. 11b. — Evolutionary synthesis model fit for source Bl and a truncated Salpeter IMF. The extinction is assumed to be Ay = 11 mag. 




Fig. 11c. — Evolutionary synthesis model fit for source C and a truncated Salpeter IMF. The extinction to the gas is assumed to be the 
same as the extinction to the stars, Ay = 15 mag. 
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Wc assume that approximately 70% of the bolometrie hi- 
minosity of B1+B2 corresponds to Bl. The best fit for 
source Bl is presented in Figure ll.b. The observed val- 
ues are well fitted at an age between 4.5 and 7Myr after 
the peak of star formation. The weak CO absorption sug- 
gests an age near the lower end of this range, ^ 5 Myr 
after the peak, or ~ 7.5 Myr from the rising half power of 
the assumed Gaussian in the star formation rate. 
Source C 

Sugai ct al. (1999) detect the 1.70 /im Hci line at about 
30% the strength of Br 10. Thus, the ionizing flux from 
this region is dominated by stars with temperatures above 
40,000 K (Vanzi et al. 1996). We have already noted the 
lack of first overtone CO bands from red giant and su- 
pergiant stars. Together, these observations require that 
the starburst in this component be extremely young, dom- 
inated by hot stars which have not yet evolved significantly 
into the red supergiant phase. However, the presence of 
significant supernova activity, as shown both by the [Fe ll] 
emission and nonthcairial radio flux, places a strong lower 
limit of ~ 3.5 Myr on the age of the starburst. 

As a first case for synthesis modeling, we consider a 
value for the extinction to the gas and to the stars of Ay = 
15 mag. From Figure 10 we get that (/dust/,/totai)K ^ 0.5. 
Again the observed CO index (value taken from Ridgway 
et al. 1994 and the conversion from Doyon et al. 1994 as 
above) and the equivalent width of Paa are corrected for 
the dust contribution to the contimmm. The best fit for 
the truncated Salpeter IMF is presented in Figure 11. c. It 
corresponds to between 3 and 6 Myr after the peak of star 
formation. Only the younger of these ages appears possi- 
ble given the evidence of extreme youth from the He i line. 
The age after the rising half power point in star formation 
is then ~ 6Myr. 

Alternately, we can correct the observables for the high- 
est possible extinction obtained from the near-infrared 
spectroscopy Ay = 37 mag (the dust emission contribu- 
tion to the X-band is less than 20% and the SNr given 
in Table 8 needs to be increased). We have to assume 
that 85% of the total bolometrie luminosity is unaccounted 
for in the observed total infrared luminosity, which argues 
against this high value of extinction. The age derived is 
even younger than for the previous fit. 

Source C 

Component C presents a dilemma. Its radio spectrum is 

nonthermal (Baan & Haschick 1990) and its flux density 
is about 50% of that of component C. However, it is very 
faint in the near infrared. Our starburst models (Figure 
11) assume a Gaussian time dependence with a FWHM 
of 5 Myr and the supernova rate leads the 2 /xm flux only 
slightly. However, as the model ages, the 2 jim flux will 
persist, so the only time to fit the properties of component 
C is during the initial period. To achieve a large ratio of 
radio to near infrared flux requires that the star formation 
episode be much shorter in duration than 5 Myr. Rather 
than reconfigure our starburst program for an instanta- 
neous burst, we use Leitherer & Heckman (1995). Their 
calculations show that, in such a burst, the supernovae 
reach their typical level abruptly starting at 3.5 Myr after 
the burst, while the red supergiants only begin to appear 
at about 4.5 Myr after the burst. Component C appears 
to be adequately explained, then, by a burst of nominally 
~ 3 X 10'^ M© with an age of 4 Myr. The excellent fit 



achieved is summarized in Table 9. 

It is often found that nuclear starbursts can only be fit- 
ted assuming vigorous star formation over 5-10 Myr, as 
we indicated in this work for component A and as studied 
in detail for M82 by Rieke et al. (1993) and Satyapal et 
al. (1997). It is interesting to discover for component C, 
a region of similar luminosity but not in a potential well, 
that a very short duration event is required. A similar sit- 
uation may hold for component C and possibly for Bl (if 
it is not the second galaxy nucleus), but they are not at 
the exact stage in their evolution that allows us to draw 
such a conclusion. 

4.6. Other Star Forming Components 

Source B2 

In addition to these four sources, we deduce that B2 is the 
site of recent star formation because of: 1.) its strong CO 
absorption (this work); 2.) its mid infrared excess (Gehrz 
et al. 1983; Keto et al. 1997); and 3.) its nonthermal 
radio emission (Gehrz et al. 1983, Zhao et al. 1997). We 
can put a lower limit of ~ 11 Myr on the age of this in- 
cident because of the very strong CO absorption feature, 
similar in depth to that of source A. As can be seen from 
Figure 11, the supernova rate for a given stellar popula- 
tion remains nearly constant from an age of ~ 6 Myr to > 
15 Myr. Scaling by the radio emission, we conclude that 
about 3 X 1O''M0 must have been converted into stars in 
B2. 

Supergiant Hil Regions 
Our Paa images reveal a large number of extremely lu- 
minous H II regions outside the components listed above. 
We have no way to correct the Paa luminosities for ex- 
tinction, so the derived values are lower limits. We have 
set a lower limit of Log (Paa luminosity) > 38.7 erg s~^, 
implying Log (Ha luminosity) > 39.7erg s~^, equivalent 
to the ionizing luminosity of 30 Doradus. We have located 
19 Hii regions in the Arp 299 system above this limit. Al- 
though we cannot model these regions in detail, we can 
use the evolutionary models to estimate roughly that each 
of them involves ~ 10^ Mq of recently formed stars, if the 
stars form according to the truncated Salpeter IMF and 
the regions are young. 

We can compare this result with the study of 30 nearby 
galaxies in the Ha line by Kennicutt, Edgar, & Hodge 
(1989). Our method of extracting ionizing line fluxes 
should give values systematically similar or perhaps lower 
than used in that work (Kennicutt, private communica- 
tion). They show that the presence of even one or two 
Hll regions at the 30 Dor level is uncommon in normal 
galaxies. Although they flnd that the incidence of super- 
giant H II regions increases for distorted and interacting 
galaxies, the extreme number of such objects in Arp 299 
is remarkable. 

Zhao et al. (1997) called attention to a region with 
bright H92a extending to the SE of component A. To a 
significant extent, this feature appears to be resolved by 
the NICMOS images into a few supergiant regions. How- 
ever, the orientation of the feature in the radio maps indi- 
cates there may be additional emission in less spectacular 
condensations. 

Massive Stellar Clusters 
Much of the past star formation in Arp 299 appears to have 
been organized into luminous clusters, similar to those 
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Fig. 12. — Histogram of the absolute H-band magnitudes of tlie faint sources detected in Arp 299 (not corrected for extinction). Also 
shown the range of H-hand magnitudes for globulaj: clusters and other IR luminous galaxies. 



found in other interacting and highly luminous infrared 
galaxies. The absolute if -band magnitudes for the clus- 
ters in Arp 299 (not corrected for extinction) range up to 
Mh = — 17mag (sec histogram in Figure 12). This distri- 
bution does not include most of the low-luminosity sources, 
since those were not detected in the NICl FllOM image. 
We can compare with the young stellar clusters found in 
the Antennae, Mh = —11 mag to Mh = —16 mag (Whit- 
more & Schweizer 1995) assuming V — H ^ 2. The derived 
sizes for the stellar clusters in the Antennae are 12 pc (for 
Ho = 75km s'^ Mpc-\ Whitmore & Schweizer 1995), 
whereas our diffraction limit at 1.1 /xm is approximately 
18 pc. In the merger galaxy NGC 7252 young clusters 
range from Mh = —16 mag to AIh = —18 mag (Whit- 
more et al. 1993). The luminosities of the stellar clus- 
ters in other infrared luminous and ultraluminous galaxies 
(Scoville et al. 1999) are also similar to those in Arp 299. 
However, it should be noted that the ULIRGs in Scoville 
et al. (1999) paper are more distant than Arp299, and 
despite the similarities in ii-band luminosities, we may be 
probing different spatial resolutions. 

The cluster luminosities in these extreme galaxies may 
exceed the limits found in more normal conditions. For 
example, the intcrmcdiatc-aged clusters in MlOO (Ryder 
& Knapen 1999) have Mh = -12 mag to Mh = -15 mag 
assuming (H — K 0.2). Since the infrared luminosities 
change only slowly (sec Figure 11), we can compare di- 
rectly to see that these values may be about 2 magnitudes 
lower than for ultraluminous galaxies. Further data are 
needed to confirm this trend. 

A crude estimate of the masses of the stellar clusters in 
Arp 299 can be obtained using the evolutionary synthe- 
sis models presented in the preceding section. Although 
at ages > 25 Myr, the luminosity will fade, it seems un- 
likely that the clusters in Arp 299 are significantly older 
than this value, given the young ages derived for the rest 
of the star formation. In Figure 13 {left panel) we present 
a color-magnitude diagram (absolute iJ-band magnitude 
versus the {H — K) color) for the stellar clusters, not cor- 



rected for extinction. The lines are the output of models 
with a truncated Salpeter IMF and masses of 1.5 x 10^ M©, 
7.5 X 10^ Mq, and 5 x 10'^ Mq. The peak of the Mh mag- 
nitude on each curve corresponds to an age of 10'' yr, the 
age at which the red supergiants make the largest con- 
tribution. The evolution is shown from 1 x 10® yr up to 
5 X 10^ yr. The star-symbols are values for the stellar clus- 
ters. The errors are 0.15 mag in the colors, and 0.1 mag 
in the absolute magnitude, and represent the uncertainties 
associated with the subtraction of the underlying galaxy 
(Section 3.2). Most of the stellar clusters lie to the right 
of the model predictions indicating that some extinction 
is present. 

From Table 5 it is clear that most of the stellar clusters 
do not show extremely red colors as it is the case of the 
bright nuclei A, Bl, B2 and CI. However the observed col- 
ors are still redder than the values predicted by the model. 
The predictions of the model show variations in the J — H 
color of less than 0.1 mag aroimd J — H = 0.7 for ages 
> 10 yr. Since the ages of the clusters are expected to be 
greater than lO'' yr, one can use the observed J—H color to 
obtain an estimate of the foreground extinction. We find 
that the extinctions are typically Ay — 1.7 — A mag, much 
lower than the extinction to the bright sources. Figure 13 
{right panel) shows how the values of Mh and {H — K) 
corrected for extinction lie closer to the model predictions, 
although they are still somewhat redder. The masses of 
these clusters are in the range of 1.5 x 10^ — 5 x 10® M©. 
For comparison the masses of the stellar clusters found in 
other infrared luminous galaxies are 10^ — 10^ M© (Scov- 
ille et al. 1999), whereas the masses of Galactic globular 
clusters arc; of the order of 10'^ Mq (van der Bergli 1995). 

Given the modest extinctions to the clusters, their deep 
CO bands (see Table 5) require little correction. Assuming 
roughly solar metallicity, these band strengths can only be 
produced by stellar populations of age ~ 10-30 Myr. 

4.7. An Extended Period of Elevated Star Formation 
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Fig. 13. — Left panel Absolute H-hand magnitude versus (H — K) for the stellar clusters, not corrected for extinction. The lines are the 
output of models using a truncated Salpeter IMF, and cluster masses of 1.5 X 10^ M0, 7.5 X 10^ Mq, and 5 x lO^ Mq. Right panel Same but 
using colors and magnitudes corrected for extinction. The error bars are indicated in the upper right corners of the plots. 



The major starforming components in Arp 299 show a 
large range of ages. Even though some uncertainties re- 
main regarding the ages and burst strengths, we find an 
age gradient from B2 (> 11 Myr) to A (~ 11 Myr) to Bl (~ 
7.5 Myr), C 5 Myr), and C (4 Myr). Compared with 
Satyapal et al. (1999), we derive an older age for source 
A, assuming that its starburst has had a relatively long 
duration (they assumed an instantaneous burst), a similar 
age for source C, and we resolve source B into two bursts 
of differing ages, one older than and the other younger 
than their value. The effective star formation rate (mass 
of the newly formed stars over the FWHM of the burst) 
is high for component A (140 M0 yr~-^). It is moderately 
high also for Bl (40 M© yr'^) and for C (22 Mq yr'^) 
(or 120 M0 yr~^ if the extinction to C is as high as 37 
magnitudes) . The masses in newly formed stars are about 
7 X 10^ Mq for component A, about 1 x 10^ Mq for Bl, 
1 X 10^ Mq for C, 1 X 10^ Mq for C, and 0.3 x 10* Mq 
for B2. 

We find a similar range for star formation in very mas- 
sive stellar clusters. We have identified 19 supergiant Hll 
regions, and have deduced total stellar masses for them 
that are similar to the masses of 21 stellar clusters with 
-15 > M^f > -17. The Hii regions must be young, 5 Myr 
or less, whereas the strong CO absorption in the clusters 
implies typical ages > 10 Myr. 

Star formation has been occurring contemporaneously 
throughout the two galaxies over the past 10-15 Myr, both 
within and outside of their nuclei. There appear to be four 
potentially differing environments for the star formation. 
Outside of the nuclei and along the spiral arms, there are 
massive stellar clusters and Hil regions. Within the nu- 
clear potential wells, young stars have formed in numbers 
one to two orders of magnitude greater than in these clus- 



ters. Although we cannot be sure of the spatial distribu- 
tion of star formation in the nuclei now, it is plausible that 
these populations will merge into a single super-cluster, 
given the high density of young stars and the strong tidal 
forces that would cause clusters to disintegrate rapidly. 
In the interaction region, components C and C appear 
to be locations where star formation has occurred on the 
scale typical of a galaxy nucleus, not an extranuclear clus- 
ter. Presumably, there is no nuclear potential well in this 
region, but the concentration of molecular gas in the in- 
teraction has resulted in an exceptionally vigorous star 
formation episode. This is also found in other interacting 
galaxies (e.g., VV114 Prayer et al. 1999; the Antennae 
Vigroux et al. 1996). Finally, there may be examples 
of extreme star formation near the galaxy nuclei where in- 
falling gas clouds became unstable and formed stars before 
falling all the way into the nuclear potential. One possible 
example is the region SE of the nucleus of IC 694 (compo- 
nent A) in the II92q! measurements of Zhao et al. (1997) 
and partially resolved into supergiant H 11 regions in our 
data. Another possible example is components Bl and 
B2. It is not obvious which should be associated with a 
galaxy nucleus. However, we have found that both involve 
a larger mass of young stars than can be attributed to a 
supergiant Hii region or massive stellar cluster, so it is 
plausible to assume one is a galaxy nucleus and the other 
a very nearby star forming region. 

5. FATE OF THE STARBURST REGIONS 

Sources Bl & B2 
The measure of stellar dynamics in B2 (Shier et al. 1996) 
suggests a mass of 6 x 10* Mq, whereas from the radio 
luminosity we have found that only about 3 x lO"" Mq is 
expected from recent star formation. Thus, this object 
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appears to satisfy expectations for a modest mass galaxy 
nucleus. On the other hand, the ^^CO emission is centered 
on Bl and has a FWHM of 260 km s'^ (Aalto et al. 1997), 
possibly indicating a larger mass there, ~ 3 x 10^ M0. It is 
possible that the CO is subject to non-gravitational forces 
at Bl, or that the stellar population at B2 is not gravita- 
tionally bound, so the situation of these two nearby con- 
centrations is not firmly determined by the observations. 

If Bl has a mass of 3 x 10^ Mq, then the Roche ra- 
dius for B2 is ^ 200 pc, about one third of the projected 
separation of the sources. In this case, B2 will be under- 
going significant tidal evaporation and is in the process of 
merging into the general stellar distribution around Bl. 
On the other hand, the dynamical mass of B2 is similar 
to the mass of Bl indicated for the recent star formation. 
We may be seeing two galaxy nuclei of comparable masses, 
one older and the other newly formed, orbiting each other 
at a projected distance 5-6 times the Roche radii. The dy- 
namical evolution will be more complex in this situation 
and should develop more slowly, but will probably have a 
similar end point. 

Sources C & C 
Sources C and C are well isolated from the galaxy nu- 
clei in the system - they are about 3 kpc from component 
A and 1.5 kpc from B, projected onto the sky. They do 
not appear to lie within any significant potential wells of 
their own, based on the lack of evidence for an underlying 
concentration of old, red stars. This case is particularly 
strong for C, given how faint it is in the near infrared. 
Both objects have low velocity dispersions in ""^^CO of 60 
- 80 km s~^ (Aalto et al. 1997). Given our estimates of 
the total mass being converted into stars (Table 8), the 
low velocity dispersion is consistent with the formation of 
a bound stellar system. The relaxation time for such a 
system can be estimated from the expression in Spitzer 
(1987) to be long, of order 10^^ yr. Thus, it would appear 
that C and C represent the formation of stellar systems 
on the scale; of dwarf galaxies. A further confirmation of 
this possibility comes from the star formation efficiency. 
The gas masses of both C and C are ~ 3 x 10^ Mq (from 
the -'^^CO measurements of Aalto et al. (1997) and our 
adjusted conversion to H2 mass, see Section 4.2). The 
mass of newly formed stars is likely to be Ix 10^ Mq for 
both, as indicated in Table 8. Thus, a significant fraction 
of the molecular material in these regions has probably 
been turned into stars, a situation that favors formation 
of bound systems. 

The efifects of tidal forces on this system can be esti- 
mated roughly from the Roche limit, which is about 500 
pc relative to the nucleus of IC 694. Therefore, the stel- 
lar systems may persist for a significant lifetime, so long 
as they do not approach more closely to the massive nu- 
cleus of IC 694 than at present. They are perhaps best 
described, given their masses and potential stability, as 
new dwarf galaxies. 

Supergiant Hll Regions & Stellar Clusters 
Although the giant stellar clusters seen in other starbursts 
have been described as young globular clusters, the distri- 
bution of these objects along spiral arms argues against 
a literal interpretation of this view. In IC 694, many of 
these objects lie a few Roche radii from the galaxy nu- 
cleus in projected distance. If the larger mass estimate 
holds for Bl in NGC 3690, a similar situation holds there. 



It is therefore plausible that the clusters that lie relatively 
close to the; massive galaxy nuclei will be subject to strong 
tidal evaporation and will fade in prominence as the mas- 
sive stars die and the less massive ones are lost into the 
surrounding galaxies. However, clusters that lie relatively 
far from nuclear mass concentrations may be reasonably 
stable. 

6. CONCLUSIONS 

We have used 775T/NICMOS broad-band and narrow- 
band images of Arp 299, along with archival HST data 
and ground-based spectroscopy and imaging to study this 
complex system. We find: 

1. ) Arp 299 is dominated by widespread and powerful star 
formation. In common with other recent studies, we find 
no evidence for an AGN anywhere in the system. 

2. ) As found for other infrared star burst and luminous and 
ultraluminous galaxies, the extinction toward the dom- 
inant infrared sources is extremely strong and complex. 
Studies using optical or near infrared measures of extinc- 
tion (e.g., Satyapal et al. 1999) will significantly under- 
estimate its importance and will deduce erroneous values 
for the critical parameters characterizing the star forming 
activity. Studies based only on UV light will face even 
more severe difficulties. 

3. ) The infrared [Fell] line emission predicts a rate 
of supernova activity in close agreement with that from 

the nonthermal radio fluxes, and both arc in satisfactory 
agreement with the predictions of evolutionary synthesis 
star burst models. 

The detailed morphology of the components of Arp 299 
has been revealed with the i/S'T'/NICMOS broad-band im- 
ages. The system is a veritable museum of interaction- 
induced star formation. It includes: 

1. ) A massive galaxy nucleus (IC 694, source A) under- 
going an extended-duration starburst and accounting for 
3 X 10"'^^ L0, about 50% of the total luminosity of the sys- 
tem. 

2. ) A strongly star forming region near this nucleus, H92a 
SE, probably indicating a region where infalling molecular 
clouds broke up into star forming regions before getting 
all the way into the nuclear potential well. 

3. ) A pair of sources (NGC 3690, Bl and B2) proba- 
bly representing a less massive galaxy nucleus and a vig- 
orous star formation episode where molecular clouds be- 
came unstable in its vicinity, together representing about 
1.5 X 10" Lq, about 27% of the total. 

4. ) Very young and vigorous star forming regions (C and 
C) near the interface between the two galaxies. They are 
likley to be forming gravitationally bound stellar systems 
that will eventually have the appearance of dwarf galaxies. 

5. ) A numerous population of supergiant Hii regions and 
stellar clusters both near the more powerful star forming 
regions and along spiral arms. They appear to represent 
young and old phases of extreme star formation outside 
the galaxy nuclei. 

We have used evolutionary synthesis modeling and other 
arguments to assign ages to these various manifestations 
of star formation, finding that they have all persisted for 
at least 10 — 15Myr. Thus, Arp 299 has been in a condi- 
tion similar to what we observe now for about this period 
of time, with all or nearly all the varieties of star forming 
environments active. 
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Arp 299 should be similar to other luminous infrared 
galaxies. The greater complexity revealed for it is primar- 
ily a result of its relative proximity (42 Mpc) and our use of 
very high resolution near infrared images from NICMOS, 
giving a physical resolution of about 20 pc. 
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